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Abstract 

We calculate the first vertex correction <5Ak to the bare vertex A Q in the 
nearly antiferromagnetic spin fluctuation Fermi liquid theory of the cuprate 
superconductor YBa2Cu307. It is calculated for k on the Fermi surface, 
and Q = (±^,±^). We find that the dimensionless ratio |<5Ak|/A , which 
parametrizes the vertex correction, is not small. It is a maximum for k = kh, 
where kh is a "hot spot" on the Fermi surface. This makes large quantitative 
corrections to the theory. 
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Most theories of the cuprate superconductors are either (i) phenomenological attempts 
to tie together many different experiments, or (ii) attempts at a "microscopic" derivation, 
beginning from models such as the Hubbard or t-J Hamiltonians. 

Of the variety of phenomenological theories, one of the most widely-discussed is that in- 
volving antiferromagnetic (AFM) spin fluctuations in a random-phase approximation (RPA) 
model of a nearly AFM Fermi liquid, in the optimally- doped regime |]-[| . For weak doping 
a Non-Linear Sigma (NLcr) model is used 0-^], but no attempt is made to deal with the 
intervening metal-insulator transition. This transition plays a central role in the more mi- 
croscopic theories, many of which do not yield a Fermi liquid state except at high T and/or 
rather large doping; instead one finds that "singular interactions" are generated and 
these inevitably give a non- Fermi liquid normal state |fFHT0|j, even at optimal doping. How- 
ever they also yield strong AFM spin fluctuations fT0| , [TT|| , as do some other microscopic 
approaches; a wide variety of phenomenological approaches involving AFM spin fluctuation 
is possible. 

To decide between competing theories it is useful to check their internal consistency. 
Several analyses have been made of the critical properties of AFM spin fluctuations JTT|-|r3||, 
but here we wish to focus specifically on the Fermi liquid model We find, in agreement 

with earlier work |[T2|Jl3[| , that once the nearly AFM form of the spin fluctuation spectrum, 



x(q,t*0 = , i c2( — 7§r 2 — —/ — C 1 ) 
i + r(q- Q) - ico/cosF 

has been assumed for q near Q = (±~,±~), and provided we assume from the outset a 
Fermi liquid model, then no singularities strong enough to destroy Fermi liquid behavior 
emerge in the fermion spectrum, provided £ and usf are finite. However we also find that 
a detailed quantitative determination of vertex corrections shows that previous estimates 
of these have been over-optimistic - they are not small. This has important quantitative 
repercussions for the theory. 

We start from a model Hamiltonian 
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^ = H e p^ PCT VV + f Yl ^k+qa^k/3(Ta/3.S(-q) (2) 



P.°" q,k,a,/3 

where 



e p = — 2t[cos(p x a) + cos(p y a)] — At' cos(p x a) cos(p y a) — fi (3) 



is the quasiparticle dispersion relation PJT41, with t = 0.25eV and t! = —0.451 The presence 



of t' allows "hot spots" on the Fermi surface Sp (FIG. [!]) which can be connected by Q; 
this leads to singular behavior when the gap in (|l|) disappears [[131 ,1511 . S(r) is the spin 
fluctuation density operator; use of (fj), with a bosonic spin fluctuation propagator ([I]) and 
a phenomenological interaction g, allows one to develop the usual "paramagnon" effective 
field theory . The first vertex correction can be calculated for general external momenta 
17]]; here we concentrate on the interaction between a Fermi surface electron and a spin 
fluctuation with q = Q (FIG. |2|) for which the lowest vertex correction is 

^ = ~jV £ f ^§X(Q' + q, «)G(k + Q' + q, u)G(k + Q" + q, Q) (4) 

where <5A k = <5A(k, e = 0; q = Q, Q = 0) and Q" = Q + Q'; the sum over Q' takes care of 
Umklapp processes. In the reduced Brillouin zone Q" = (0,0). Note that the overall sign 
of the graph in FIG. |2] is negative, ie., eqtn. (^) is negative; the bare vertex A Q = § in eqtn. 
(0) is positive. 

We concentrate on <5A k , with k e Sp (here Sp is the Fermi surface) because it gives an 
indication of the size of vertex corrections. This is because in general |5A(k, e; q, fl) | exceeds 
|5Ak| (in fact it diverges along a surface in k-space, for given values of q, Q and e); thus one 
cannot argue, even after integrating over one or more of its arguments, that |<5A(k, e; q, Q)\ 
will lead to corrections smaller than one would get from just using <5A k (our argument here 
parallels Migdal's fl8||). More generally one finds that if Q <C qvp, A , where A is the spin 
gap, then 5A k is a good approximation to 5A(k, e; q, fl). Writing (f|) as 

6A ^__9^yY^f d2(1 dU X"(Q' + <1,U) ( /k+Q'+q , I" /k+Q'+q 



A 4 q, 7 (27r) 2 7r e k+ q/ +q — e k+ Q» +q e k+ Q/ +q — u e k+ Q/ +q + to 

/k+Q"+q l-/k+Q'+ Q N 2 XQ (USF^n- (r\ 



Ck+Q"+q - U ^k+Q'+q + U 12tC 3 UJsF V 
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where g 2 = |g 2 and X k is dimensionless; we define g to correspond directly with the coupling 
constant g used in Monthoux and Pines ||. At zero temperature one has 

/2-7T r27T I 
dq x j dq y - — (Sgn(ei)Gi(q, ei) - Sgn(e 2 )Gi(q, e 2 ) - G 2 (q, ei) + G 2 (q, e 2 )] 
~2?r J-2tt e± — C2 

(6) 

where q x = q x a, etc.; ei = e k+Q ' +q /yU and e 2 = e k+Q » +q //i = e k+q /yU. d and G 2 are defined 
by 



(S) 



withX=[l + (i) 2 ^ + g2)](c 5F /H). 

X k can be investigated both analytically and numerically Here we calculate it for two 
different points in k-space, both on the Fermi surface (see FIG. |1|); ki makes a 45° angle 
with k x , and kh is a "hot spot" wave vector. 

In order to obtain a value for <5A k , we need values for the spin fluctuation energy losfi 
the correlation length £, the susceptibility XQi the coupling constant g, and the chemical 
potential \x (which is determined by the electron filling factor n). There are different values 
reported in the references fl3 Ui"9|j2Cfl . We have evaluated <5A k for k = k 1; kh, in two ways, viz. 
(a) by assuming various published values for the different parameters, and (b) by making 
the simple assumption that one is very close to an AFM instability, and then, in the spirit 
of RPA, imposing the condition |(^/2)n( 1 )(Q, 0)| = 1 where n«(Q , uS) is the electron-hole 
bubble. Numerical calculation gives |II«(Q,0)| = 2.6(eV)" 1 and thereby g = 0.67 eV, 
assuming n=0.75 and the band structure in (|3|). Since u)sf/\iA ^ 1> this value of g should 
be a very good guess, within a naive RPA scheme. 

The results are summarized in Table |I[ We use two different values for ujsf ; these are 
the two different values quoted by Monthoux and Pines et al. ||. We also use two different 
values for the coupling constant g, quoted from MPI and MPII (see ref. again). We use 



values of £ = 2.5a, \Q — 80 states/eV (from |J) and n=0.75 , appropriate to YBa2Cus07 
(again quoted from [|J). This value of n corresponds, with the band structure in (0), to a 
value of \y\ ~ 1.46* = 0.365 eV. 

We see that even the values for the vertex correction calculated from the simple RPA 
model (b) are not small; as in the standard discussion of Migdal's theorem, the importance of 
vertex corrections appears in the ratio |<5Ak|/A . If one takes values of g from the literature 
||, this ratio is quite unreasonably large (as large as 2.43 for the hot spots in the model 
used by MPI). Thus vertex corrections are clearly very important. The values we quote 
for |o"Ak|/A D are considerably larger than previous estimates [p, 22|,23| . The reason for this 



difference with previous work can be tracked back to the factor 1^, which is impossible to 
guess from purely dimensional arguments. In fact if we drop the factor Xy_ from 5A^, we get 
an order of magnitude estimate for SA-^ given by 

|MkL -O^T^TT 1 « 1 (9) 



A 47T 3 |/i| \fl 

which is broadly in agreement with previous estimates (see eg. Millis |22|) ; xq ; ^sf and g 
must be redefined to conform with the parametrizations in this paper). 

In fact however 1^ is surprisingly large, and also shows a significant variation around the 
Fermi surface, with a maximum at the hot spots, and a minimum at intermediate wave- vector 
like ki. We should emphasize here that analytic calculations of SA^ have to be approximated 
rather carefully in order to give reasonable agreement with the numerical results in Table 
|. Approximations such as those of Hertz et al. J24j (see also p3j), which try to separate 
off a rapidly- varying (in q-space) contribution from %"(q, a;), give quantitatively incorrect 
results (including a completely unphysical ln[(k — kh)a] divergence as one approaches the 
hot spot). A more detailed discussion of the behavior of is given in ref. |T7f . 

One might suppose that X k is large simply because of the band structure (ie., because 
of van Hove singularities, or the hot spots). If this were true one could argue that the 
quasiparticle weight ought to be renormalised down near these singular points in the Brillouin 
zone, and that this would considerably reduce the vertex correction. In fact however we find 



this is not the case; this can be checked analytically by suppressing the regions immediately 
around the hot spots in the integral for 1^, or by simply redoing the numerical calculation for 
a slightly different band structure. We find that suppressing the hot spots entirely, reduces 
the vertex correction by a factor which is everywhere less than 2 (and which differs very 
little from unity when k is far from a hot spot). Thus we do not believe that incorporating 
self- energy corrections near the hot spots would significantly reduce the vertex correction. 

We re-emphasize here that these results do not necessarily invalidate the internal con- 
sistency of the Fermi liquid starting point, in this theory. However they do show that the 
theory cannot be trusted quantitatively, at least in the usual RPA form. As is well known 
the RPA is not a "conserving approximation", and for spin fluctuation theories this makes 
it unreliable (cf. ref. fl~6fl , especially section 3). It is useful to compare the case of nearly 



ferromagnetic 3 He liquid, where vertex corrections are also quite large, and where use of 
the paramagnon model yields values for m*/m which are off by a large factor [^1]]. Thus 
if we use melting curve Landau parameters, ~ 0.75 and Ff ~ 15, we infer a value for 
the Stoner factor S ~ 24 which yields m*/m = ^\nS ~ 15 , in the paramagnon model. 
This is roughly 2.5 times the correct value of ~ 6 (note that the first vertex correction is 
S A/ 1 ~ InS ~ 3 in this model), and no amount of self-consistent summing of diagrams can 
cure this numerical problem. 

Similar problems can clearly occur in the present AFM spin fluctuation model. We believe 
this is the main reason for the difficulty one encounters in the MP models, in determining 
a value for g that (a) gives the correct superconducting T c , and (b) is consistent with the 
observed spin susceptibility. 

To check the structure at higher order, we have also estimated the contributions from 
the graphs containing 2 spin fluctuation lines (there are actually 7 distinct graphs at this 
level) , and found that some of them are also large for the values of g used above |L7| . Thus, 
just as for the case of nearly ferromagnetic 3 He, we see no reason to believe, for the values of 
the parameters given in the table, that performing infinite graphical sums will lead to results 
which are numerically more reliable, even if they do converge to some smaller renormalised 



vertex-there will always be other diagrams with large values, which will in general give 
uncontrolled contributions. 

It is interesting to compare these results with some other recent investigations. In the 



weak-coupling limit, Chubukov |13| has calculated the leading vertex corrections to g, con- 



centrating on the gapless case; in the case where is a gap, he finds that the renormalised 



coupling is also small (for realistic values of the gap). On the other hand Schrieffer [25 
has argued that a correct formulation of the theory, even in the weak-coupling limit, must 
take account of the short-range local antiferromagnetic order even in the normal state-if this 
done, he finds that a weak-coupling calculation shows very strong suppression of the vertex 
when one is close to the antiferromagnetic transition. This theory seems rather interesting- 
note that a related calculation by Vilk and Tremblay finds that the existence of such 
a short-range antiferromagnetic order will cause a breakdown of the Fermi-liquid starting 
point itself! Thus the question of what is the correct theory itself is rather confused, even 
in the weak-coupling regime. It is certainly not clear how any of these arguments will work 
in the regime discussed in this paper, when g is not small enough to control the magnitude 
of the vertex corrections. 

It is of course crucial that these higher-order corrections also be included in any version 



of this theory that tries to reconcile different experiments - as emphasized by Pines pO 
the justification of the theory stands or falls on its ability to do this quantitatively. It is 
possible that such a programme might succeed if one can show that the actual parameters 
g, usf, and xq are such that |<5A k |/A is considerably less than one (ie., if one is genuinely 
in the weak-coupling regime). This would also be true of versions of the theory in which 
usf depends on g, whilst the spin gap becomes an independent parameter [|13|]; or of the 
theory of Schrieffer cited above [25|| . On the other hand if |5Ak|/A > 0(1), we see no hope 
that such a scheme could succeed quantitatively (in., eg., the calculation of T c ), since the 
vertex corrections become large. 

We would like to thank I.Affleck, V.Barzykin, D.Bonn, and W.Hardy for useful conver- 
sations, as well as A. Chubukov for correspondence. This work was supported by NSERC in 



Canada, and by the CIAR. 
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TABLES 





9 (eV) 


u>SF (meV) 




5A kl /A 




<5A kh /A 


MPI 


1.36 


7.7 


78.6 


-1.81 


105.6 


-2.43 


MPII 


0.64 


14 


49.6 


-0.46 


73.4 


-0.68 


"RPA" 


0.67 


7.7 


78.6 


-0.44 


105.6 


-0.59 


14 


49.6 


-0.50 


73.4 


-0.74 



TABLE I. Calculated values of the vertex correction <5Ak for two different wave-vectors ki and 
kh on the Fermi surface (columns 4 and 6 in the table). The papers MPI and MPII (ref. g) 
give different values for g, and different values for u>sf- From the values for these two models one 
calculates Zk in equation (|6|), and thence 5Ak/A . The third model is the naive "RPA" model 
described in the text, for which g is determined; we have calculated 2k and <5Ak/A for two values 
of ujsf given in MPI and MPII respectively. 
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FIGURES 




FIG. 1. Fermi surface in the first Brillouin zone, with the value for t and t' given in the text; 
and we assume n=0.75. Calculations are presented here for the wave vectors k\ and k^. 
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FIG. 2. First correction 5A^ to the bare vertex, for an incoming fermion with momentum k 
and energy (relative to the Fermi energy), interacting with a fluctuation of wave- vector Q and 
zero energy. 
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